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We report on a preliminary study of the volume dependence of the gluon propagator. The propaga- 
tor is computed using different lattice volumes, its extrapolation to infinite volume is investigated 
with particular attention to the its IR behaviour. Our data shows a mild but measurable depen- 
dence with the volume. Unfortunately, we are not able yet to clarify its behaviour close to zero 
momentum. 
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1. Introduction and motivation 

The infrared (IR) properties of the gluon propagator were studied in [jl], ^ with the following 
large assymetric lattices 16 3 x 128 and 16 3 x 256. For these lattices, although the propagator for 
pure temporal momenta are within errors, the equivalent spatial momenta show a deviation from 
the temporal ones for small momenta. Clearly, this is a finite volume effect. So the question, to 
which we will try to give a first and preliminary answer, is how reliable are the figures reported 
in [fl], |2p? Or what is the meaning of these results? In order to try to answer these questions, we 
make a first study of the volume dependence of the gluon propagator, which will also allow an 
extrapolation to infinite volume. 

2. Lattice setup 

In order to study the volume dependence we generated 1 SU(3) pure gauge, Wilson action 
configurations for the lattices reported in table [I]. The table also describes the combined over- 
relaxed+heat bath Monte Carlo sweeps, the number of combined sweeps for thermalization, the 
separation of combined sweeps used for each simulation and the total number of configurations for 
each lattice. 



Lattice 


Update 


therm. 


Sep. 


Conf. 


8 3 x 256 


70VR+4HB 


1500 


1000 


80 


10 3 x 256 


70VR+4HB 


1500 


1000 


80 


12 3 x 256 


70VR+4HB 


1500 


1000 


80 


14 3 x 256 


70VR+4HB 


3000 


1000 


47 


16 3 x 256 


70VR+4HB 


3000 


1500 


155 



Table 1: Lattices setup used in the study of the volume dependence. 

For the definition of the gluon propagator, notation and the gauge fixing method see [gp. 

3. Time versus Spatial Momenta 

Our first check was to look at how the time and the different spatial momenta evolved with the 
lattice volume. The gluon propagator, for the smallest and the largest lattices, for different types 
of spatial momenta is given in figure [T]. The data shows that the propagator becomes smaller when 
the number of zero components increases. The effect is larger for smaller momenta. For the largest 
lattice, the effect is not so strong but is, nevertheless, clearly visible. 

The comparison between temporal and spatial momenta with only one nonzero component 
shows a similar behaviour. The temporal data is always below the spatial data - see ph. 

The gluon dressing function for temporal momenta and spatial momenta with only one nonzero 
component show a strong dependence on the volume, see figure ||, specially in the intermediate 
momenta region (~ 1 GeV). Note that while the temporal data seems to increase with the volume, 

'All configurations were generated with MILC code http : / /physics . indiana . edu/ " sg/milc . html. 
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Figure 1: Dressing function for different spatial momenta, for the smaller (left) and largest (right) lattices. 



the spatial data decreases with the volume. Certainly, this is the case for almost all the momenta 
region but not for the lowest momenta - see figure [|. Indeed, a zoom to the IR region shows the 
opposite behaviour, i.e. for small momenta the propagator decreases with the volume. Therefore, 
if in the IR the gluon dressing function is given by (q 2 ) lK , a decreasing of the dressing function 
means that K should increase with the volume. Then, the values measured with any finite lattice 
can be seen, at least, as a lower bound in K. 
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Figure 2: Gluon dressing functions 



4. IR region 

In [Q| it was seen that the IR region is well described either by a pole or a cut behaviour 2 . 
We fitted both expressions to the different lattices. As seen in figure ||, there is good aggreement 
between the two expressions. Excluding the smallest lattice, the K (certainly not A) seems to be 

2 See [fj] for details. 
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stable against a change in volume. For the largest volume K = 0.5138^22 an d 0.5122^}j for the cut 
and pole formulas, respectively. In this sense, the lattice data supports a vanishing zero momentum 
gluon propagator. 



0.48 
0.46 



Figure 3: Z cut and Z„ i e fitted parameters as function of lattice volume. 



Let us now look at the behaviour of the propagator as a function of the volume. For the 
four lowest momenta values, the volume dependence can be seen in figure ^. To extrapolate the 
zero momentum gluon propagator we tried the formulas reported in table 2. The results are a 
bit disapointing in the sense that all functional formulas produce good fits. From the table, one 
sees that the extrapolation is compatible with all types of zero momentum behaviour, from zero to 
infinite. Nevertheless, the results of table 2 will help to reduce the number of fitting formulas to be 
considered in the extrapolation of the full propagator. 
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Figure 4: Propagator as a function of the inverse volume. 



287/4 



Finite volume effects in the gluon propagator 



Orlando Oliveira 





X 2 /d.o.f. 


Fitted Parameters 


a + bx 


0.60 


a = 72.06, b = 3.0 x 10 6 


a + bx + cx 2 


0.24 


a = 69.52, b = 5.1 x 10 6 , c = -2.5 x 10 11 


a + bx c 


0.10 


a = 53.96,6= 1702, c = 0.32 


ax b 


0.12 


a = 325.4, 6 = 0.11 


{a + bx)x c ' 


0.13 


a = 2 x 10~ 5 , 6 = 418, c = -0.87 



Table 2: Zero momenta propagator extrapolation, where x = 1 /V . 



5. Infinite volume extrapolation 

As a guide to the infinite volume extrapolation we take the results of the zero momentum gluon 
propagator. Since in the last fitting formula a ~ 10~ 5 , we will only consider the following fitting 
functions 

a + bx + cx 2 , a + bx c , (5.1) 

and will fit only the timelike momenta. Moreover, each timelike momenta will be fitted separatly. 
From the two fitting formulas, the second one does not provide a smooth propagator. For certain 
momenta the fits do not converge or produce an infinite propagator. Therefore, from now on we 
will consider only the fits assuming a quadratic behaviour. Indeed, for the quadratic function, the 
extrapolated propagator looks rather smooth (see figure ||) and only less than 15% of the fits have 
zX 2 /d.o.f. >2. 
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Figure 5: Extrapolated gluon propagator. 



If one tries to repeat the study performed in [p], ||] it cames that the only fits which have 
acceptable % 2 jd.o.f . are the fits to the IR region. Probably, this is due to the inclusion of the data 
from the smallest lattice which looks different from all the other ones. The results of the IR fits 
are reported in table ||. Again, the lattice although pointing towards a vanishing zero momentum, it 
does not provide a definite answer. Looking at the k values, the fit to the analytical DSE solution is 
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compatible with 0.5. However, it is interesting that if one includes corrections to the above solution, 
the result agrees within errors with the figures estimated from the largest lattice, using the same 
procedure. 





Extrapolated 

K A x 2 /d.o.f. 


16 3 x 256 
k A X 2 /d.o.f. 


(<?) 


0.4993(7) 0.37 


0.4858(2) 0.40 


(? 2 r(i + ^ 2 ) 


0.5129(10) 0.07 


0.5070(50) 0.44 


/ 2 \ 2 K 

W+A*) 


0.5198(2) 438(1) 1.29 


0.5090(20) 409(4) 0.71 


(q l Y K 


0.5167(1) 439(1) 1.36 


0.5077(17) 409(4) 0.69 



Table 3: IR fits. The errors on the extrapolated propagator are clearly underestimated. 



6. Conclusions 

We have performed a first study of the volume dependence of the gluon propagator. Although 
the lattice data seems to favour K > 0.5, we are not able yet to provide a clear answer concerning 
the behaviour of the gluon propagator at zero momentum. In what concerns the IR region, our data 
shows a mild but measurable dependence with the volume. This is not in aggreement with a recent 
similar study of the DSE equations on a tours [[|, ||] . 
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